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Abstract 
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that found in quantum optics. The chaoticity parameter is strongly momentum dependent. 
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1 Introduction 



In current ultrarelativistic heavy ion collision experiments at the SPS (e.g., Pb + Pb at Ef, eam = 160 
AGeV) secondary particles are formed at high number densities in rapidity space In future 
experiments at RHIC and the LHC one expects to obtain even higher multiplicities on the order 
of a few thousand particles per unit rapidity. If local thermal (but not chemical) equilibrium 
is established and the number densities are sufficiently large, the pions may accumulate in their 
ground state and a Bose condensate may be formed^]. A specific scenario for the formation of a Bose 
condensate, namely, the decay of short-lived resonances, was discussed in Ref. || where conditions 
necessary for the formation of a Bose condensate in a heavy ion collision were investigated. In Ref. 
H it was found that if a pionic Bose condensate is formed at any stage of the collision, it can be 
expected to survive until pions decouple from the dense matter, and thus it can affect the spectra 
and correlations of final state pions. 

In the present paper we investigate the influence of such a condensate on the single inclusive cross 
section and on the second order correlation function of identically charged pions (Bose-Einstein 
correlations BEC) in hadronic reactions for expanding sources. 

The proposal to use BEC for the detection of condensates was made a long time ago ||. In that 
reference the influence of the condensate on the form and the intercept of the second and third 
order correlation function was studied, with particular emphasis on the possible role of the strong 
final state interactions, which, contrary to more phenomenological approaches used afterwards, 
were studied using a quantum statistical Landau-Ginsburg type method. 

It was found in || that the interaction influences the shape of the correlation function but not 

the value of its maximum (intercept). On the other hand the amount of condensate influences 

both the intercept and the form of the correlation function. In particular in analogy with quantum 

optics, a typical two component structure in the correlation function was found: while for purely 
1 Recently a different type of pion condensate, the disordered chiral condensate, has received a lot of attention in 
the literature (for a review of recent results, see Ref. fy). It has been argued || that such a condensate would lead 
to the creation of squeezed, i.e., two particle coherent states. The possible effects of such squeezed coherent states on 
Bose-Einstein correlations have been studied in some detail [Q. 

In what follows, however, we shall restrict ourselves to the case of a conventional Bose condensate which is formed 
when the pion chemical potential becomes equal to the pion mass and to the associated one-particle coherent states. 
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chaotic fields the second order correlation function has only one term which depends on the rapidity 
difference, for a superposition of coherent and chaotic fields there are two such terms which differ 
in a well defined manner and which depend both on the same correlator (cf. eq,(|23|) below). No 
dependence on the total rapidity of the pair appears in this approach because of the assumed boost 
invariance and of the fact that in the one dimensional treatment (in rapidity space) used in || no 
allowance for expansion was made. 

The present paper uses a hydrodynamical approach which differs significantly from the previous 
approach in that it considers explicitly the expansion of the source and thus the correlation between 
position and momenta of produced particles. Such an approach is certainly more adequate for heavy 
ion reactions and appears desirable at this stage particularly because of the possibility mentioned 
above that a condensate might be produced after the expansion of the system during the freeze- 
out process. Furthermore, the momentum dependence of the chaoticity can easily be taken into 
account. Last but not least, in the present paper we exploit explicitly the possibility offered by 
quantum statistics to treat simultaneously single inclusive cross sections and BEC and to correlate 
the possible effects of a condensate in these observables, which is of considerable interest for future 
comparisons with experiment. 

The remainder of this paper is organized as follows. In section 2 general expressions are derived for 
the single particle spectrum and two-particle correlation function of pions emitted from a hydrody- 
namically expanding fluid with a superfluid component. In sections 3 and 4 the formalism is applied 
to the case of a spherical expansion and a longitudinal scaling expansion, respectively. Numerical 
results for both models, obtained for specific parametrizations of the freeze-out hypersurface, are 
presented in section 5. Finally, the main results of the paper are summarized in section 6. 



2 Basic formalism 

Applying the current formalism ||, ||, |l(]] we write the amplitude for the emission of a pion of 
four-momentum k = (E, k) as a superposition of contributions from the thermal excitations and 
from the condensate, 

J(k) = J th {k) + J co {k) , (1) 
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with 



*' e lkxi j\ h (k) , 



(2) 




(3) 



n n 



where the labels "th" and "co" indicate the thermal and the condensate component of the distribu- 
tion. The indices / and n label source elements centered at space-time point x\ and x n , respectively, 
and ex.Y>(i4>\ h ) and exp(i^°) are the corresponding fluctuating phase factors. 

For the thermal component, two phases that characterize the emission from two different source 
elements, (pf 1 and cfift with I / I', are taken to be uncorrelated. Likewise, we assume that there 
exists no correlation between the phases related to emission from the thermal component and from 
the condensate, 4>\ h and 4> c ° . If, and to what degree, the phase factors for emission from the 
condensate are correlated for two different source elements depends on the details of the formation 
and evolution of the condensate. There are two extreme cases: 

(a) The phases § c ° are completely correlated, i.e., the differences — § c °, do not fluctuate. In the 

of case an expanding source this would imply that the condensate must have been formed at 
an early stage when the hadronic matter was concentrated in a space-time volume sufficiently 
small for this kind of global phase coherence to be established. This corresponds to the 
presence of an expanding superfluid with a globally coherent phase. 

(b) The phases 4> c ° are completely uncorrelated. This is what one would expect if the condensate 

is formed at a late stage when the process of condensate formation (caused, e.g., by the decay 
of short-lived resonances) occurs independently in different fluid cells. 

In the following we shall consider only the case (b), i.e., we assume that the Bose condensate is 
not phase correlated in different fluid cells. This situation appears as more realistic than case (a) 
although the formation of a global superfluid as envisaged in (a) is in itself also of interest and 
might deserve further study. 

The single and double inclusive momentum distribution are 



E 



d 3 N 



= (J*(k)J{k)) 



(4) 
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and 



E\E 2 - 



d 6 N 



{r(k 1 )j*(k 2 )j(ki)j(k 1 )) , 



where the averaging is performed over the space-time positions xi, x n and the phases (frf 1 , . 
It is useful to introduce the following correlators: 



(5) 



D™{k l M)^(Y.( J n(k l )TJ™(k 2 )), 

G th (k u k 2 )^(j2\jt h (h) 2 Jf h (k 2 



(6) 
(7) 

G C °(k U k 2 )^(J2 \ J n( k l)\ 2 \Jn{k 2 )\ 2 ) . (9) 

The particle spectrum may then be expressed as the sum of a thermal and a condensate contribution, 

d 3 N 



d 3 N ^d 3 N 

E—r— = E 



with 



and 



d 3 k d?k 
d 3 N 



E 



E 



d?k 
d 3 N 



th 



d?k 



th + E d 3 k 



D th (k,k) 



D co (k,k), 



(10) 

(11) 
(12) 



and the two-particle inclusive distributions may be written as 



(J + (fe 1 )J*(fe 2 )J(fe 2 )J(fc 1 )) 



D^ik^k^+D^i^M) 
-G^ik^k^-G^ih,^). 



(13) 



The two-particle correlation function then takes the form 

(J*(fc 1 )J + (A; 2 )J(A: 2 )J(A; 1 )) 



C 2 (ki, k 2 



<J*(fcl)J(fcl))(J*(fc 2 )J(fc2)> 



1 + 



D m {k l , k 2 ) + D co (k u k 2 ) - G m (h,k 2 ) - G co (h,k 2 ) 



{D th {ki, fci) + D co (ki, k!)){D th {k 2 , k 2 ) + D co (k 2 , k 2 )) 



(14) 



We first consider the thermal correlators D th (ki,k 2 ) and G th (k\,k 2 ) to further evaluate these 
expressions. Let L be the number of source elements that contribute to particles of momenta 
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k\ and k 2 . In the limit of large L the term G th {k\,k 2 ) oc L can be neglected compared to the 
term \D th (ki, k 2 )\ 2 oc L 2 in eq. In the case of hydrodynamics, the sum over source elements 

translates into an integral over the freeze-out hypersurface. Introducing the average and the relative 
four-momentum of the pair, K = \{k\ + k 2 ) and q = k\ — k 2 , respectively, one has 

D th (k 1 ,k 2 ) = [ d 4 x g th (x,K) e iqx , (15) 



with the thermal source function for pions [11] 



g th (x,k)~ !h 



(2vr) 3 7e 



k da(x') 8 4l {x — x' 


) 


exp 


k u(x') — Hir(xf) 


- 1 


Tf(x') 



(16) 



In eq. (16) da(x), u(x) and Tf(x) are the volume element of the freeze-out hypersurface E, the 
4- velocity of the fluid and the freeze-out temperature at space-time point x, respectively. The factor 
g n denotes the degeneracy factor of the pions while fi n (x) is the pionic chemical potential at space- 
time point x. For simplicity, in the following sections we shall assume that /jL n (x) = /% = const., i.e., 
that a pionic Bose condensate is formed at each space-time point on the freeze-out hypersurface. 

Having discussed the terms D th (ki, k 2 ) and G th (ki,k 2 ) which are due to the thermal part, we 
now proceed to consider the corresponding expressions related to the condensate, D co (ki,k 2 ) and 
G co (ki, k 2 ). In what follows each fluid element will be treated as a macroscopic system in so far as 
it will be assumed that the condensate in each fluid cell is identified with the lowest momentum 
state in the rest frame of that fluid cell. This is in line with the conventional meaning of the concept 
of a condensate which refers to a phase of matter. That is to say, for a fluid cell centered around 
xi moving with velocity ui we have 

\Jf°{k)\ 2 oc 5 3 (k- m7T ui), (17) 

i.e., as emphasized in Ref. ||, particles emitted from the condensate move with the collective 
velocity of the fluid. In the models of expanding sources which will be discussed in Sections 3 and 
4 there is a one-to-one relation between the position x\ of the source element and its velocity u\. 
As far as particles emitted from the condensate are concerned, eq. ( |i~7|) then implies that for each 
momentum k only one single source element contributes to the spectrum. Note that nevertheless 
this contribution will be comparable in magnitude to the contributions of the thermal excitations 
since the condensate constitutes a macroscopically occupied quantum state. 
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From eqs. ©, (|) and (0) it follows that 



D co (k 1 ,k 2 ) = G co (k 1 ,k 2 ) = forfeit. 

The contribution of the condensate to the single inclusive spectrum can then be written as 

d?N 



d?k 



D co (k,k) = I d A xg C0 (x,k) 



with H 



g co (x,k) = / (la(x')u(x') n co (x') E 5 3 (k — m^u) 5 4 (x — x') . 
Jt, 



(18) 



(19) 



(20) 



We define the fraction of thermally produced particles of four- momentum k, i.e., the momentum 
dependent chaoticity p(k), as 

D th (k,k) 



E 



p(k) 



d 3 k 



tii 



E 



d 3 N D th (k,k) + D c °{k,k) 
IFk 



(21) 



and the normalized correlator of two thermal currents 

D th (k 1 ,k 2 ) 



22) 



[D^{kiM)D th {k 2 M)] 1/2 ' 

Substituting the expressions (|15|) , ( ^6[ ) and ( |i~9| ) - ( p2| ) into (Q) we find for the two-particle Bose 
Einstein correlation function 



c 2 (h,k 2 ) = i + 2^/ P (k 1 ) P (k 2 )(i- P (k 1 ))(i-p(k 2 )) e{\h-k 2 \) d th (h,k 2 ) 

+ p(k 1 )p(k 2 )\d th (k l ,k 2 )\ 2 , (23) 



where 



Q{x) 



i.e., one has 



1 : x < 
: x > 

th i 



C 2 (k 1 ,k 2 ) = l+p(k 1 )p(k 2 )\d th (k 1 ,k 2 )\ 2 for k x + h 



Note that the "true" intercept - obtained by evaluating (^3|) at k\ = k 2 - is given by 

C 2 (k, k) = l + 2p(k) -p\k) , 
whereas the value obtained by extrapolating (p5|) to k\ = k 2 is 



(24) 



(25) 



(26) 



a 



extrapol. 



(k,k) =l+ P \k) 



(27) 



The reason for the difference between (|26| ) and (|27|) is that in the idealized case considered here 
D co (ki,k 2 ) and G co (A?i,fe 2 ) vanish for k\ + k 2 (cf. Eq. (||)), while 
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G co (k,k) = D co (k,k) 



(28) 



where I labels the single source elements which contributes to the emission of pions of momentum 
k from the condensate (cf. eqs. (@||0))- 



3 Spherically expanding source 

Let us assume a space-like freeze-out hypersurface parametrized as t = t(r) where t = x° is the time 
coordinate and r = \x\ the radial coordinate. We also assume a radial velocity field u = u(r) e r 
with non-vanishing gradient, i.e., (du/dr) ^ 0. The pionic freeze-out temperature is taken to 
be Tf = m n = const. The volume element of the 3-dimensional freeze-out hypersurface and the 
4- velocity of the spherically expanding relativistic fluid are (u = 0, 1, 2, 3) 

da v {x) = (l, ^e r ) r 2 dr sin 6 d9 dcp , (29) 



dr 

u u {x) = (u°(r),u(r)e r ) , u°(r) = y/l + [u(r)] 2 . (30) 
3.1 Single inclusive momentum distributions 

With the 4-momentum k v = (E , k) = (E , k ■ et) and e r ■ et = cosO = z, the thermal part of the 
single inclusive momentum distribution is given through 



d 3 k 



^ do v k v 



th ( 27r ) 3 exp [(u u k u - fJhr)/Tf] - 1 

j2^Io ^^I-i dZ exp [(Eu°(r) - kzu(r) - n w ) /T f ] - 1 ' (31) 

In eq. (ETI) E = \J m'^ + k' 2 and R± is the radial extension of the pion source 



The coherent part of the single inclusive momentum distribution is 
d 3 N 



E 



d?k 



n co / do y u u E8 (k — m^u) 



E 



f - 1 - f n f ( dt 

n co / f 2 dr I sin# d6 I dcp I u° (r) — u(r) — 
Jo io j o V or 

x i? 5 3 (kek — m n u(r)e r ) 

^£(£-^fM)(f) 2 (f(ro)) _1 : 

. k ^> TTT'Tv^max 



max 



(32) 



where ro = r(u = ^-) and ti ma x is the maximal value of the velocity at freeze-out. The total 
momentum distribution is obtained by inserting eqs. (H) and §D into eq. ©. 



3.2 Bose-Einstein correlation functions 



For illustration we restrict ourselves to the central momentum region. To be specific, we consider 
the case K = -^(k\ + fo) = 0, i.e., k\ = —k~2- The quantity that remains to be calculated in 
order to construct the BEC function is the thermal correlator D {k\, fe). With E\ = E2 one has 
q° = Ei — E2 = 0. Using the notation q = k\ — k2 we obtain q v x v = —q - x = —qrcosO and 

D th (k 1 ,k 2 " 9w 1 v 



K=o (2vr) 3 h exp [{K v u v - ^)/T f ] - 1 

9* ,if n -n,nf R± 2r EkC^ COS t 



s'mOdO / r dr 



(2vr)3y Y J J exp[(E K u°(r)- ^)/T f ]-l 

9 * ^rdr u J E « S ™ iqr K q , , • (33) 



(2vr) 2 7 exp [(^u°(r) - ^)/T f ] - 1 " 

In section 5 explicit parametrizations for t(r) and u{r) will be used to evaluate numerically the 
above expressions. 



4 Longitudinally expanding source 

Here we assume a space-like freeze-out hypersurface parametrized by the radius-dependent longitu- 
dinal proper-time r(r). For the longitudinal component of the velocity we take the scaling ansatz 
vu = xu/t. To be specific, the 4- volume element of the freeze-out hypersurface and the 4- velocity 
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of the relativistic fluid are written as 

dr ^ 

da v = ( cosh j) , sinhr? , — e r ) r(r) dri r dr dip , (34) 

or 

u u = (U(r) cosh rj , U(r)smh. rj , u_i_(r)e r ) , U(r) = \Jl + [«j_(r)] 2 , (35) 

where uj_(r) and 77 are the transverse (radial) component of the fluid 4- velocity and the space-time 
rapidity of the fluid, respectively. 



4.1 Single inclusive momentum distributions 



Since our ansatz is invariant under boosts in the longitudinal (beam) direction we may without loss 
of generality restrict ourselves to the case of rapidity y = 0. With the 4-momentum k v = (E , k) = 
(m ± cosh y , m± sinhy , k^e^) and e r ■ = cos ip the thermal part of the single inclusive momentum 
distribution is given by 



E 



d 3 N 



d 3 k 
9- 



th 



(Jtt 



da v k v - 



v = o 



(2vr) s 



(2vr) 3 7s exp [{u v k v - ^)/T f ] - 1. 

rR x r-2-K 

dr] / r(r)rdr / dtp [m± cosh r\ — k± cos ip (dr/dr)] 
Jo Jo 

1 



exp [(mj_ cosh 7/ U(r) — k± cos ip u±(r) — fi w )/Tf] — 1 



(36) 



In eq. fl36|) R± is the radial extension of the pion source and m± = \/m% + k\. 



The coherent part of the single inclusive momentum distribution is 

d 3 N 



E 



d 3 k 



n. 



CO 



da v u v E8 {k — m^u) 



y = 



ii,,, j dr] / r(r)rdr / d(f> I U(r) — itj_(r)— I 



oo JO JO 
•3/ 



<9r 7 



x mj_ 5 (k — m^u) 



y = 



: > rn^Umax 



(37) 



where ro = r(uj_ = ^-) and u maa; is the maximal value of the transverse velocity field at freeze-out. 
The total momentum distribution is obtained by inserting ( |36[ ) and (37) into eq. (10). 
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4.2 Bose-Einstein correlation functions 



As in the preceeding section we restrict ourselves to the case K = -(^i + ^2) = ^, i.e., k\ = 
—ki- Since we have assumed the source to be invariant under boosts in longitudinal direction, the 
generalization to average pair momenta K± =0,Ku ^ is straightforward. We shall consider the 
correlation function in the transverse direction C2(q±,Ay = 0), i.e., the case K» = fcjii = yi = 



and Ek = E{ = m_\_K = ym^ + kf ± (i = 1, 2). Then for construction of the BEC function the only 
remaining term to specify is D(ki, A^)- 

With q° = Ei — E2 = we obtain q v x v = —q- x = —q±rcosip (since q\i = 0) and 



D(k u k 2 ] 



k=o (27r) 3 Jy. exp \{K p u v — fj, n )/Tf] — 1 

Q-k f +oa , f 2n f R± , T(r)m ±K cosh 7] e- ig±rcos ^ 
- (^loo dr, J ^Jo rdr exp [(E K U{r) - ^)/T f ] - 1 
= _9 ]L _ f + °° d r R ± rdr T(r)m ±K cosh 77 J (q L r) ^ 

(27r) 2 J -00 Jo exp [(m ±K U(r) cosh?? - fi^/Tf] — 1 ' 

with the Bessel function Jo. 

In the next section explicit parametrizations for r(r) and u±(r) will be given and the results for 



eqs. (|36|), (|37|) and (|38|) will be discussed after their numerical evaluation. 



5 Parametrizations and Results 

For the sake of illustration of the above results, we shall consider the case of a specific heavy-ion 
reaction, namely S + S at 200 AGeV. Application of the relativistic (3+l)-dimensional hydrody- 
namic code HYLANDER [12, [l3j yielded a successful description jnj of rapidity and transverse 
momentum spectra of mesons and baryons and a quantitatively correct prediction of BEC of pions 

MM- 

In order to get reasonable descriptions for t(r), u(r), r(r) and u±(r) we have chosen the following 
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parametrizations and parameters from the hydrodynamical solution mentioned above: 





< 




r < R± 



with Toc/fm =ao/fmc = 4.5 



(39) 







r>R ± 



u(r) = u±(r) = 



u r 



r<R± 



with uq fm/c = 0.13 . 



(40) 



< 







r > R ± 



The pionic chemical potential was taken to be ix n = 0.139 GeV, since we consider the formation of 
a pionic Bose-condensate. The transverse radius has been fixed to R± = 3.8 fm. 

The remaining parameter is the condensate number density n co which we write in fractions of the 
thermal number density 



In Fig. 1 the results of the numerical evaluations of the single inclusive momentum distributions 
EdN/d s k, the momentum-dependent chaoticities p and the Bose-Einstein correlation functions C2 
are shown for the spherically and for the longitudinally expanding source. The quantity n co here 
takes four different values ranging from 0.0 nth to 1.0 nth- 
It can be seen that the results for the spherical and the longitudinal expansion scenario look quite 
similar although there are some quantitative differences according to the momentum dependence. 
The Bose-condensate affects the plotted functions only over a limited momentum range. This is 
due to the fact that there exists a maximum velocity u max = u(R±) = u±(R±) = 0.494 c that 
enters in all condensate contributions. Thus the maximum momentum where the Bose-condensate 
contributes to the single inclusive momentum distributions is h max — f^TvUmax — 

68.7 MeV/c and 

the maximum momentum difference over which the BEC can have a reduction is twice this value, 
namely q ma x = 

137.4 MeV/c. 

Of course, theoretically the largest value u max could assume is the speed of light c = 1. There- 
fore, there exists an absolute maximum value for the momenta of particles which originate from 
a Bose-condensate: in case of a pionic Bose-condensate it would be k max = m^c = 139.6 MeV/c. 
Bose-Einstein correlation functions would be then affected over a momentum difference range with 




(41) 
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Qmax = 279.2 MeV/c, which would result in a disappearance of the sharp peak in the BEC for 
S + S, because this peak would be shifted to unobservably large momentum differences. It should 
be mentioned that from such a structure in the double inclusive signal one could in principle de- 
termine the maximum flow velocity of the fluid. 

The single inclusive momentum distributions are quite sensitive to the presence of a Bose-conden- 
sate; for the BEC this effect is even more pronounced. The presence of a Bose-condensate of about 

only 1% ( 1.6. /"v — ^co 

= 0.01 n t h) results in a decrease of the intercept of the shown two-particle 
correlation functions to 1.83 for the spherically expanding source and 1.86 for the longitudinally 
expanding source. In quantum optics as well as in the approach used in || the sensitivity of the in- 
tercept on the amount of coherence is much weaker. Thus to get a decrease of the ("true") intercept 
from 2 to 1.8, a value of k = 0.5 is necessary. Furthermore, due to a limited, value of Qmax 

a part of 

the tail of the two-particle correlation functions will not be affected by the pionic bose-condensate 
and a peak will therefore come into existence. To what extent such peaks due to an existing pionic 
Bose-condensate can be observed in experimental data depends on several factors: 

(i) the size of the source, which determines the inverse width of the correlation functions; 

(ii) the maximum of accessible velocity at freeze-out of the pionic source which determines the 
position of the peak; 

(iii) deformation of the correlation function due to further contributions from resonance decay; 

(iv) the averaging over acceptance regions in the experiment; 

(v) the width of the momentum distribution in the bosonic ground state. 



A study of the infiuenece of these factors on the effects reported here is in preparation [16] 



6 Summary 

We have shown that the formation of a pionic Bose condensate can influence single inclusive spectra 
as well as Bose-Einstein correlation functions quite strongly. Among other things we find an 
enhancement of single inclusive momentum spectra at low momenta, and a reduction of the two- 
particle correlation function on a limited momentum range resulting in a bumpy structure of the 
BEC. The structure depends on the maximum velocity of the fluid. Within this treatment an 
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explicitly momentum-dependent chaoticity is presented for the first time due to a maximum freeze- 
out velocity rather than to geometrical effects of the source. 

Whether these effects will survive in more realistic treatments has yet to be proven. However the 
results obtained so far confirm the interest and need to take into account coherence in theoretical 
approaches to BEC. 

This work was supported by the University of California and the Deutsche Forschungsgemeinschaft 
(DFG). B.R. Schlei acknowledges a DFG postdoctoral fellowship and R.W. is indebted to A. Capella 
for the hospitality extended to him at the LPTHE, Univ. Paris-Sud. 
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Figure Captions 



Fig. 1 Single inclusive spectra, momentum-dependent chaoticities and BE correlation functions 
for a spherically and a longitudinally expanding source, respectively. The different line styles 
correspond to different condensate densities n co compared to the thermal number densities 
nth- The left column of the figure shows the results for the spherically expanding source, 
whereas the right column shows the results for the longitudinally expanding source. 
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